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Adsorption of acetylene on Si(lOO) is studied from first 
principles. We find that, among a number of possible adsorp- 
tion configurations, the lowest-energy structure is a "bridge" 
configuration, where the C2H2 molecule is bonded to two Si 
atoms. Instead, "pedestal" configurations, recently proposed 
as the lowest-energy structures, are found to be much higher 
in energy and, therefore, can represent only metastable ad- 
sorption sites. We have calculated the surface formation en- 
ergies for two different saturation coverages, namely 0.5 and 
f monolayer, both observed in experiments. We find that al- 
though, in general, the full monolayer coverage is favored, a 
narrow range of temperatures exists in which the 0.5 mono- 
layer coverage is the most stable one, where the acetylene 
molecules are adsorbed in a 2 x 2 structure. This result dis- 
agrees with the conclusions of a recent study and represents a 
possible explanation of apparently controversial experimental 
findings. The crucial role played by the use of a gradient- 
corrected density functional is discussed. Finally, we study 
thermal decomposition of acetylene adsorbed on Si(fOO) by 
means of finite-temperature Molecular Dynamics, and we ob- 
serve an unexpected behavior of dehydrogenated acetylene 
molecules. 



I. INTRODUCTION 

Adsorption of unsaturated hydrocarbon molecules on 
Si surfaces is of great importance to investigate the 
initial-stage growth of SiC films on Si substrates. In 
particular, there is an increasing interest in the study 
of adsorption of two basic hydrocarbon units, namely 
acetylene, XC2II2) and ethylene (C2H4), on Si(lOO) 
surfaces.ErO Acetylene, due to its reactive triple bond, 
has a high probability of interacting with the surface 
Si dangling bonds, undergoing electronic rehybridization 
and sticking to the surface. Adsorption of acetylene 
on Si(lOO) represents a promising carbon source for ho- 
moepitaxial growth of cubic SiC because the processing 
temperature can be reduced well below 1000 °C, thus 
making it compatible with the present silicon technolo- 
gies. In fact, at temperatures as high as 600 °C, acetylene 
already decomposes, leading to desorption of hydrogen 
and the formation of SiC clusters on the Si(lOO) surface. 
Furthermore, the diffusion of a small concentration of 
carbon atoms into the Si substrate may form a substoi- 
chiometric alloy with interesting IR optical properties. 

For the adsorption of acetylene on Si(lOO), as in the 
case of ethylene, most of the experimental and theo- 
retical studies support the "di-cr" model, in which the 
acetylene molecule is adsorbed across a first-layer Si 
dimer in a bridge (B) site, with the Si dimers pre- 
served (i.e. not cleaved) upon adsorption. Evidence 
favoring thej-.|B-| ajdsorption site comes from different 
experimentsJda'aEl based on high-resolution electron- 
energy-loss spectroscopy, low-energy electron diffrac- 
tion. X-ray photoemission spectroscopy, near-edge X- 
ray-adsorption fine structure, scanning tunnelling mi- 
croscopy, phatpelectcon diffraction, and from theoretical 

calculations .Ereliaia 

Recently, hcevever, this scenario has been challenged 
by Xu et al.fB who studied the adsorption of ethylene 
and acetylene on the Si(lOO) surface by using a photo- 
electron holographic imaging technique. The most sur- 
prising result of this study was the observation of a new 
preferred adsorption site for acetylene on Si(lOO), where 
the molecule is bonded to four silicon atoms, in a pedestal 
configuration (P) between two adjacent Si dimers. This 
so-called "tetra-a" model seems to bax;onsistent also with 
high-resolution photoemission data. El 

By analyzing their diffraction images, Xu et al. infer 
a C-C distance in adsorbed acetylene of 1.2 A. Even 
taking into account the indeterminacy associated to this 
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estimate (which, according to the authors, could be as 
large as 0.2 A), this represents another surprising result. 
In fact this value is comparable with the experimental 
bond length of the triple C-C bond in the free acetylene 
molecule (1.21 A), whereas, for the sp^^-hybridized acety- 
lene adsorbed in the P site, one would expect a much 
longer bond, close to the single C-C bond of, say, C2H6 
(1.54 A). 

Another somewhat controversial issue concerns the 
C2H2 equilibrium coverage which is most likely to be 
realized on the Si surface. Based on experimental mea- 
surements, both 0.5 monolayer (ML) and 1 ML cover- 
ages have been proposed as saturation coverages real- 
ized, in equilibrium conditions, on the Si(lOO) surface (1 
ML means that as many C2H2 molecules are adsorbed 
as Si surface dimers). In particular, STM experiments 
performed at room temperature indicated that the C2H2 
molecules form a local p{2 x 2) or c(4 x 2) structure with 
a saturation coverage of 0.5 ML. Detailed analysis of 
the STM images reveals lack of adsorption on nearest- 
neighbor sites, thus showing that the C2H2 molecules 
adsorb on alternate dimers: the two observed surface pe- 
riodicities depend on whether the C2H2 molecules ad- 
sorbed on adjacent rows are in phase or out of ahase. 
On the other hand kinetic uptake measurements ,0 per- 
formed at lower temperatures (105 K), provide evidence 
for 1 ML coverage, i.e. each Si dimer site adsorbs a sin- 
gle C2H2 molecule. These different exp^mental obser- 
vations have been tentatively explained^ by assuming 
that the early stages of acetylene chemisorption are ki- 
netically rather than thermodynamically driven, i.e. ad- 
sorption proceeds via precursor states,Ll which although 
not truly thcrmodynamically stable states, are separated 
by sizeable potential barrier from the stable adsorption 
configurations. 

The process of thermal desorption of acetylene on 
Si(lOO) is also of interest. It has been foundtl that, as 
the surface temperature is raised up to 870 K, the acety- 
lene molecules dissociate, leading to the formation of SiH 
groups and surface carbon. However, the details of the 
dissociation paths and the early stages of surface car- 
bonization are not known. Simulations performed us- 
ing ab-initio finite-temperature Molecular Dynamics can 
provide valuable information to elucidate such complex 
phenomena. 

In order to clarify the above issues, we have performed 
a full ab initio study of acetylene adsorption on Si(lOO), 
by using the Car-Parrinello approach.ES 

The paper is organized as follows. In Section II we 
describe the computational method. We report our re- 
sults for the adsorption geometries, with particular em- 
phasis on the energetic and structural characterization 
of the B and P configurations, in Section IIIA. In Sec- 
tion IIIB we address the problem of the theoretical de- 
termination of the C2H2 equilibrium coverage. Our cal- 
culations of the surface formation energies for different 
model structures at 0.5 and 1 ML coverage provide evi- 
dence that the saturation coverage depends on tempera- 



ture, largely favoring the 1 ML coverage, but indicating 
the existence, for a particular adsorption configuration, 
of a narrow window in temperature (which we find only 
qualitatively in agreement with the experiment), where 
the 0.5 ML coverage may become favored. We offer this 
prediction, which disagrees with the conclusions of Ref. 
p^ , as a semi-quantitative explanation for the different 
experimental observations quoted above. 

Finally, in Section IV we present the results of our 
simulations of thermal dissociation of C2H2 on Si(lOO); 
this is a very important topic being related with the pro- 
cess of growth of SiC films obtained from decomposition 
of acetylene chemisorbed on the Si(lOO) surface. Our 
results are encouraging and show some interesting be- 
havior, including a tendency of the acetylene molecules 
to decompose into intact C dimers, which diffuse below 
the Si surface layer and bind themselves to form more 
complex structures. 

II. METHOD 

Total-energy calculations and Molecular Dynamics 
(MD) simulations hp^e been carried out within the Car- 
Parrinello approaches in the framework of the density 
functional theory, both in the local spin density approx- 
imation (LSDA) and using gradient corrections in the 
BLYP implementation.^ Gradient-corrected functionals 
have been adopted in the most recent theoretical stud- 
ies of adsorption of organic molecules on Si(lOO) because 
they are typically more accurate than the LSDA func- 
tional in. -describing chemical processes on the Si(lOO) 
surfaceEJa. The calculations have been carried out con- 
sidering the F-point only of the Brillouin-zpne (BZ), and 
using norm-conserving pseudopotentialsJ£3 with s and p 
nonlocality for C and Si. Wavefunctions were expanded 
in plane waves with an energy cutoff of 50 Ry. We have 
explicitly checked that, at this value of the energy cutoff, 
the structural and binding properties of our system are 
well converged. 

In order to test the above approximations, we have 
calculated the structural and vibrational properties of 
the free acetylene molecule (see Table | ) and find an 
excellent agreement with the exppimental estimates and 
previous theoretical calculations.ES 

The Si(lOO) surface is modeled with a periodically re- 
peated slab of 5 Si layers and a vacuum region of 7 
A (tests have been also carried out with a vacuum region 
of 10 A, without any significant change in the results). 
A monolayer of hydrogen atoms is used to saturate the 
dangling bonds on the lower surface of the slab. We have 
used a p(4 x 4) surface supercell in the (2 x 1) reconstruc- 
tion, corresponding to 16 Si atoms/layer (some calcula- 
tions have been performed using instead a (4 x 2) recon- 
struction, see Section IIIB). With such a relatively large 
supercell we expect that a sampling of the BZ limited 
to the F point is adequate. A sensitive test for the ade- 
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quacy of the fc-point sampling is provided by the study of 
the clean Si(lOO) surface: as discussed below, the correct 
surface structure is obtained with our supercell. 

Structural relaxations of the ionic coordinates are per- 
formed using tbj2 method of direct inversion in the itera- 
tive subspace.L^ During ionic relaxations and MD simu- 
lations the lowest Si layer and the saturation hydrogens 
are kept fixed. We verified that, by starting with the un- 
reconstructed, clean Si(lOO) surface, the structural opti- 
mization procedure correctly produces asymmetric sur- 
face dimers, with a dimer bond length and buckling an- 
gle in good agreemen|L-|With previous, highly converged 
ab initio calculations.L3 Acetylene molecules are added 
on top of the slab and the system is then fully relaxed 
towards the minimum energy configuration. To better 
explore the complex potential energy surface of this sys- 
tem, in most of the cases the optimization procedure was 
repeated using a simulated-annealing strategy and also 
starting from different initial configurations. The same 
optimization procedure has been already successfully ap^ 
plied to the study of adsorption of benzene on Si(100).til 

^^o^S nr-i 
the many possible adsorption configurationg3il3 for acety- 
lene on Si(lOO) we have focused our study on the B and 
P structures (those involved in the present controversy 
about the lowest-energy site), although calculations have 
been performed for other configurations as well. Adsorp- 
tion of acetylene between two adjacent dimer rows is not 
considered in our calculations, since it is known ta,p:sult 
in structures having much lower binding energies .Hlij 

III. RESULTS 
A. Lowest-energy adsorption sites 

Four possible configurations of acetylene adsorbed on 
Si(lOO) are shown in Fig. |l|; the energetics and the struc- 
tural properties are summarized in Table ||. When com- 
paring the total energies of different adsorption config- 
urations, we consider a representative 0.5 ML coverage 
of acetylene on Si(lOO), corresponding to two Si dimers 
per C2H2 molecule. This is, the coverage realized in the 
experiments by Xu et aLf3 and the one which is likely 
to be prevailing at room temperature. Note that calcula- 
tions performed at lower coverage (for instance at 0.125 
ML coverage, corresponding, in our supercell, to a single 
acetylene molecule adsorbed on the slab) do not alter the 
conclusions contained in Table |l[ as far as the energetic 
ordering of the structures is concerned. 

We confirm that both the P and B structure are pos- 
sible adsorption configurations, and that the Si dimers 
involved in bonding with acetylene are preserved (in line 
with the most recent experimental and theoretical works, 
see Refs. |6|,p^ and further references quoted therein), al- 
though the typical buckling of the (2x1) reconstruction 



is removed and the Si dimers become symmetric (this lat- 
ter result is also in agreement with previous calculations, 
see, for instance, Ref. |l6|). According to our data the P 
structure, however, lies at much higher energy than the B 
one. Interestingly, a lower-energy pedestal configuration 
(P ') is obtained by rotating the P structure by 90° with 
respect to the Si surface. In this configuration, the four 
C-Si bonds are shorter than in the P one. However, P ' 
is also found to be less favored than B. If the B structure 
is rotated by 90°, in such a way to bridge two Si atoms 
of two adjacent dimers, one obtains another possible ad- 
sorption configuration, B ', slightly higher in energy than 
B. Our results compare favorably with those of the recent 
paper by Sorescu and Jordan ,Ej where however the ad- 
sorption energies for the B and B' structures at 0.5 ML 
coverage are not reported. We also find that, at higher 
coverage, the B ' structure becomes favored with respect 
to the B one (the energy difference is ^ 0.14 and 0.19 
eV/molecule, using the LSDA and BLYP functional, re- 
spectively), again in_agreement with the calculations of 
Sorescu and JordanO 

As can be seen in Table use of BLYP gradient cor- 
rections instead of LSDA makes bond lengths about 1-2 
% longer; moreover the P ' configuration is much less 
energetically favored, and the B and B ' structures are 
almost degenerate. Note however that the energetic or- 
dering of the four different structures remains unchanged. 

Our results contradict the findings of Xu et alEl 
that the P site is the stable adsorption site, and con- 
firm instead the conclusions reported in the previous 
literature. For the B configuration, which is cer- 
tainly the favored one at coverages not higher than 0.5 
ML, our structural parameters are in excellent agree- 
mcnt_both with the experimental estimates of Matsui 
et alB, (C-C bond length=1.36 ± 0.04 A) and Terborg 
et alB (C-C bond length=1.36 ± 0.19 A, C-Si bond 
length=1.83 ± 0.04 .A-), a n d .-with the most recent the- 
oretical calculations.ll3"ll30"ll3 

In the B structure our optimized C-C bond length 
(1.35-1.36 A) is close to that of ethylene (1.33 A), sug- 
gesting sp^ rehybridization of theX^-C bond, sp^ bonding 
is supported by Imamura et al.,c3 who reported charge- 
density plots where the characteristic 7r-bonds are ob- 
served, and by Matsui et al.p who observe tt* features in 
near edge X-ray adsorption fine structure spectra. Our 
computed vibrational frequency of the C-C stretching 
mode is 1440 cm~^, which compares favorably both with 
the experimental value of Huang et aln (estimated to 
be in the range 1450-1500 cm~^ ) and, with the previous 
ab-initio estimate by Imamura et a/£j (1479 cm^^). 

Note that, as expected, the C-C bond length in the 
P and P ' structures is significantly larger than in the B 
one, and is not consi SLtent with the value, 1.2 ± 0.2 A, 
observed by Xu et al!Ei. This leads to the conclusion 
that, either the experimental indeterminacy associated 
to the C-C bond length is largely underestimated, or the 
structure observed by Xu et al. is not a pedestal one. 

For the B structure, again at 0.5 ML coverage, our 
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computed adsorption energy (78 kcal/mol using LSDA, 
65 kcal/mol with BLYP) is-jij-ffie«|i|agreement with pre- 
vious theoretical estimateijil30't3 (ranging from 55 to 
77 kcal/mol). The substantial energy difference between 
the bridge (B, B') structures and the pedestal (P, P') 
ones strongly supports the conclusion that the bridge ad- 
sorption sites are the stable configuration for acetylene on 
Si(lOO), while the pedestal structures can represent only 
metastable adsorption sites. 

Other possible adsorption structures exist, besides 
those shown in Fig. |l] (see also Ref. 19). For instance, we 
found a "diagonal-bridge" (DB) configuration in which 
the acetylene molecule forms a cross-bridge structure 
with two Si atoms belonging to two neighbour Si dimers 
of the same row. This configuration can occur as one 
considers reconstructions of the Si surface involving al- 
ternating buckled Si dimers, such as the p{2 x 2) and the 
c(4 X 2). However, according both_tp our calculations 
and to those of Sorescu and Jordanjlj the DB structure 
is even energetically less favored than the P ' one. 



B. Saturation Coverage 

Having assessed the quality of our approach by com- 
paring our structural and energetical data with previous 
theoretical and experimental studies, we are now in the 
position to deal with the issue of saturation coverage, 
where a high accuracy is essential to get meaningful in- 
formation. |-| 

STM measurementsQ of C2H2 on Si(lOO), performed 
at room temperature, indicate a saturation coverage 
of 0.5 ML (one acetylene molecule adsorbed every 
two Si surface, dimers). However, other experimental 
measurementaJ suggest instead a saturation coverage of 
1.0 ML. We show in the following that, although the 1.0 
ML coverage seems to be generally preferred, in prin- 
ciple both coverages can be realized depending on the 
temperature of C2H2 deposition and on the adsorption 
configuration. 

The 0.5 ML coverage can be realized in a number of 
ways, even assuming that each of the acetylene molecule 
is adsorbed in the lowest energy structure, i.e. the B 
structure (see Fig. ^). In particular, as shown in the 
following, the relative energy ordering is influenced by 
the structure of the Si dimers which are not saturated 
by C2H2 molecules. These may be characterized by a 
different relative buckling orientation. Here we adopt the 
same notation used in Ref. |l^ to identify these different 
structures. The full monolayer (1 ML) coverage, on the 
contrary, is compatible only with a single configuration if 
the acetylene molecules are adsorbed in the B structure. 

Using first-principles calculations Tanida et alx3 pre- 
dict, for the most stable 0.5 ML coverage structure, one 
(structure (a) in Fig. |^) where alternated dimer rows are 
fully covered with C2H2 molecules. However such struc- 
ture is not compatible with the STM observations of Ref. 



^, performed at room temperature. According to these 
observations, at 0.5 ML coverage the acetylene molecules 
are instead adsorbed on alternate dimers in the same 
rows, and on adjacent rows they can be in phase or out 
of phase. Note that instead the (c)-(f) configurations, 
which are characterized by (2 x 2) or (4 x 2) periodicity, 
are all compatible with the STM observations of Ref. E . 

Using the BLYP functional we confirm (see Table |l) 
the result of Ref. |l^ that the (a) structure is at the low- 
est energy. However, at variance with Ref. |l|, we find 
that the (a), (c), and (e) structures are almost degener- 
ate in energy. This, discrepancy could be due to the fact 
that Tanida et altB use a smaller cell than ours and pos- 
sibly to the different gradient-corrected functional they 
adopt. Note that in this case the effect of the gradi- 
ent corrections is not a simple quantitative correction to 
the LSDA results; in fact using the BLYP functional the 
energetic ordering of the different configurations is com- 
pletely changed (see Table |l|) . 

At thermal equilibrium the optimal surface structure 
can be determined by computing the grand canonical 
potential. t3 We show in Fig. ^, for the structures (a) 
and (e) of Fig. |[ the behavior of AO, that is the dif- 
ference between the grand canonical potential of a given 
structups and that of the clean Si(lOO) surface. Afl is 
definedliS as 



An = AE- Nr. 



(1) 



where AE is the difference between the total energy of 
the considered structure and that of the clean Si(lOO) 
surface (using the lowest-energy surface reconstruction 
which is compatible with the given configuration, that is 
(4 X 2) for configuration (a) and (2x1) for configuration 
(e), respectively), and iV „ „ is the number of the C2H2 
molecules; following Ref. [isfthe zero of the chemical po- 
tential A'caHa assumed to be 1/3 of the energy of the 
isolated benzene (CgHg) molecule in vacuum. Calcula- 
tions were performed using the BLYP gradient-corrected 
functional. As can be seen, for structure (a) the preferred 
saturation coverage appears to b<p-4-0 ML, in agreement 
with the findings of Tanida et aL;ll3 and the same is true 
for the structures (b), (c), (d), and (f) (the correspond- 
ing Ar2(/i^^jj^ ) curves are not shown because they closely 
resemble those of structure (a)). However, in the configu- 
ration (e), there is instead a range of values of the chem- 
ical potential for which the 0.5 ML saturation coverage 
becomes favored. Note that, considering the adsorption 
structure B' in place of B, the 1.0 ML saturation cover- 
age would be always favored, since this structure is the 
preferred one at high coverage (see previous discussion). 

Configuration (e), with a 2 x 2 periodicity, is compat- 
ible with STM measurements at room temperature and, 
according to our calculations (see Table ||), is also one 
of the lowest-energy configurations at zero temperature, 
being almost degenerate with configuration (a); since it 
is the only structure for which the 0.5 ML coverage could 
be favored instead of the 1.0 ML coverage, we have tried 
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to estimate in which range of temperatures this particu- 
lar configuration is stable. To this aim we have followed 
the procedure described in Ref. |3[ The surface forma- 
tion energy, whose minimum determines the equilibrium 
state of the surface as a function of composition, is de- 
fined as fl/A, A being the surface area and 



(2) 



is the grand canonical potential (we have omitted a PV 
term which, for the pressures considered here, is negligi- 
ble) . F is the free energy of the system, iV^ is the number 
of particles of type i, and fii denote the chemical poten- 
tials of the various components. If one assumes that the 
Si(lOO) surface is in equilibrium with bulk Si (this is the 
natural reservoir for Si atoms, due to the presence of 
steps, terraces and other surface defects), then for the 
bare Si(lOO) surface one has: 



ribarc = i;[Si(100)] - 7Vsii^bulk[Si] 



(3) 



Our calculations show that, in the temperature range 
from 410 to 460 K, the adsorption configuration (e) at 
0.5 ML coverage is energetically favored with respect to 
that corresponding to 1 ML coverage. Since this temper- 
ature interval is not far from room temperature (given 
the approximations involved in our derivation the accu- 
racy of our temperature estimate is of course limited) we 
conclude that configuration (e) could be realized as a sta- 
ble adsorption configuration at room temperature. This 
result, which disagrees with the conclusions of Tanida et 
al.iB (according to Ref. |l^ the 0.5 ML coverage is never 
favored), is in cmalitativc agreement with the experimen- 
tal observationd that full monolayer coverage is realized 
at T = 105 K, while the 0.5 coverage is observed in STM 
experiments,^ performed at room temperature. Of course 
this does not rule jOuJ. the possibility that the existence 
of precursor statesQO could play a crucial role in desta- 
bilizing the surface species for coverages greater than 0.5 
ML; however it offers an alternative/complementary ex- 
planation of the experimental findings. 



where the first term is the total energy of the slab and the 
second term is the energy per atom of bulk Si multiplied 
by the total number of Si atoms of the cell. In the above 
expression the temperature dependence of the free energy 
of the slab and that of the bulk reservoir have been ne- 
glected, as usually done in this kind of calculationsrj The 
formation energy of the acetylene-covered surface can be 
estimated, following the same reasoning, by computing 

n{T) = £;[Si(100):C2H2] - A^sii?buik[Si] - iV^^^, A^c.h, (r) + 

(4) 

where -Fkds is the free-energy contribution due to the vi- 
brations of the adsorbed C2H2 molecules. In order to es- 
timate the temperature dependence of /^cHa assume, 
quite reasonably, that the surface is in equilibrium with a 
reservoir of acetylene molecules in the gas pkase. There- 
fore, using elementary statistical mechanics,0 

/^c^H, (T) = £^c,H, + kT\n{PVQ/kT) - fcTlnZ^t - fcTlnZvib 

(5) 

Here -EC2H2 is the total energy of the isolated acetylene 
molecule, P is the pressure, and Vq — {h^ / 277 mkT)^^^ 
is the quantum volume. The rotational and vibrational 
contribution to the chemical potential are calculated us- 
ing the experimental values of the vibrational frequencies 
for the free acetylene molecule.El We estimate the pres- 
sure P from realistic growth conditions of C2H2 films 
on Si(lOO): experimentallyaa the saturation coverage is 
realized under an acetylene exposure of ~ 2 L, which cor- 
responds to P ~ 2 X 10~^ Torr. The vibrational contri- 
bution Fads(T) = EJ^^»/2 + kT\n{l - expi-HcuJkT))] 
is estimated from the vibrational frequencies calculated 
in Ref. ^ (the frequencies of the hindered modes,EJ T^,, 
Tj, and R^, which are not included among the mode fre- 
quencies reported in Ref. 



10, have been estimated here 



by using a "frozen-phonon" approach). 



C. Thermal decomposition of C2H2 on Si(lOO) 

As a result of the high desorption activation energy 
of C2H2 on Si(lOO), adsorbed acetylene molecules tend 
to be retained on the surface up to relatively high tem- 
peratures. According to the, EELS and LEED mea- 
surements of Nishijima et al.p upon heating the C2H2- 
EsB^Tt^, Si(lOO) surface from 80 up to 750 K, the C-H 
bond scission occurs, with a gradual recombination of 
the dissociated H atoms onto the Si(lOO) surface; the C- 
C bond scission takes place by further heating the sample 
at 870 K, while by heating up to 930 K the surface hy- 
drogen is probably completely removed by H2 desorption; 
similar qualitative results haveJpeen obtained in theoret- 
ical calculationS|-by Zhou et aLEl 

The proposedB decomposition path is therefore: C2H2 
C2ilx (a;=0,l) — * C ^ /3-SiC; the activation harrier 
for the dissociation of a single H atom is estimatedQ to be 
no less than 0.8 eV, while that for the the dissociation of 
the remaining C2 species is about 1.3 eV. This indicates 
that acetylene favors dehydrogenation to C2Ha; {x—0,1) 
species first; then the remaining C2 dimers further de- 
compose into C species to form /3-SiC films. However 
it is not at all clearDO if, after acetylene dehydrogena- 
tion, all the H atoms desorb as a gas of H2 molecules or 
if instead some of them actually form Si-H bonds. Also 
the weakening of the Si dimer bonds, as a consequence of 
heating the sample, could be relevant: according to Ima- 
mura et al.,t2l although the dimerized structure is more 
stable than the dimer-cleaved one, the latter structure 
may be taken instantaneously just after the adsorption, 
and the Si-Si bond-breaking and rebonding oscillations 
could occur at relatively high temperatures. Finally, in 
view of the SiC-growth characterization, it would be very 
interesting to clarify how, after acetylene dissociation. 
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carbon atoms diffuse into the bulk of the Si crystal. 

We have performed Car-Parrinello MD simulations at 
finite temperature in order to investigate the dissociation 
of acetylene on Si(lOO). The ionic degrees of freedom 
have been integrated using a time step of 5 a.u. (~ 0.12 
fs) and typical MD runs were 3-4 ps long. We start with a 
configuration where the C2H2 molecules are adsorbed on 
the Si(100)(2 x 1) surface in the B structure at 0.5 ML 
coverage. The initial geometry is optimized using the 
method described in Section II. We gradually heat this 
structure by simply rescaling the ionic velocities. Obvi- 
ously we are forced to use heating rates much higher than 
in actual experiments in order to let the system undergo 
some kind of transformations; this however would lead 
to melting of the whole sample prior to C2II2 desorp- 
tion/decomposition. In order to overcome this unphysi- 
cal effect, we have adopted a simple z-dependent heating 
scheme, in which the temperature decreases as one moves 
towards the interior of the slab. In this way the ad- 
sorbed acetylene molecules and the uppermost Si atoms 
are kept at higher temperatures than the remaining ions, 
thus avoiding melting of the whole slab. Although both 
in the simulations and in actual experiments the details 
of the heating procedure may significantly affect the fi- 
nal system configuration, we believe that our scheme is 
able to reproduce, at least at a qualitative level, the basic 
features of the desorption process. 

The main results of our MD simulations can be sum- 
marized as follows (given the heating procedure followed, 
our "simulation temperature" is not particularly mean- 
ingful for comparison with the experimental one): 1) by 
heating the sample the first relevant event is the dissocia- 
tion of some C-H bonds (typically one H atom is detached 
from every molecule); 2) then C-Si bonds are continously 
broken and reformed leading to diffusion of the acetylene 
molecules on the Si surface; 3) by further increase in the 
temperature, some C2II2 molecules loose their H atoms, 
remaining on the surface as C-C dimers; 4) these dimers 
remain intact even at very high temperatures; 5) the dis- 
sociated H atoms prefer to form new bonds with Si atoms 
rather than forming II2 molecules, although in our simu- 
lation this tendency could be enhanced by the finite size 
of the vacuum region; 6) the C-C dimers eventually form 
complex sub-surface structures, just below the first Si 
layer, such as a 5-atom ring composed by one Si atom 
and 4 C atoms. 

Some of the above results are consistent with experi- 
mental observations. In fact, in the early stages of an- 
nealing, a strong SiJL. stretching mode is visible in the 
HREELS spectrum,aa thus supporting the formation of 
Si-H bonds, while the peak due to C-C stretching broad- 
ens and develops a long tail; this latter result could in- 
dicate that breaking of C-C dimers could be only partial 
or/and much slower than H dissociation. The surpris- 
ing stability of C-C dimers has been confirmed by using 
different heating rates and also adopting the LSDA func- 
tional instead of BLYP. This remarkable effect is in con- 
trast with previously proposed desorption mechanisms 



which just assume a simple C-C dimer thermal break- 
ing (see for instance Ref. |l]). Note that the same result 
has been obtained by Cicero et alS3 from simulations 
in which they studied the heating of carbon dimers ad- 
sorbed on Si(lOO). 

This however poses the problem of how SiC (whose 
formation is believed to require breaking of C-C dimers) 
is actually produced. A possible explanation could be the 
presence, in real experiments, of a gas of-silanes (such 
as Si2ll6) which could play a crucial roleEJ in the SiC 
growth: in fact adsorbed Si atoms can react with the 
C2II2 molecules, and C-C dimer breaking could actually 
occur before acetylene adsorption. 



IV. CONCLUSIONS 

Adsorption of acetylene on Si(lOO) has been investi- 
gated using ab initio simulations. The structural and en- 
ergetic properties of possible adsorption configurations 
have been studied. Although the P structure, recently 
suggested by Xu et al. as the lowest-energy configura- 
tion, is found to be stable, our results show that the 
bridge structures (B or B', depending on the coverage) 
are energetically favored, in agreement with most of the 
previous studies. The results reported by Xu et al. could 
be influenced by the presence of defects, such as steps, 
which could make the actual Si surface used in their ex- 
periment substantially different from the ideal Si(lOO) 
surface assumed (see Fig. 2 of Ref. ^0|) in the interpreta- 
tion of the experimental data. 

We have then addressed the issue of the saturation 
coverage and its dependence on the temperature, and we 
find that, although the full monolayer coverage is gen- 
erally preferred, there is a temperature window, not far 
from room temperature, in which a particular adsorption 
structure at 0.5 ML coverage is preferred. This could ex- 
plain the results of STM observations performed at room 
temperature. The effects of including gradient correc- 
tions have been also discussed: our calculations show 
that these effects are particularly relevant for the deter- 
mination of the optimal adsorption configuration at 0.5 
ML coverage, where even fine structural details, such as 
the relative orientation of the buckled Si-dimers, may be 
important for the correct determination of the surface 
formation energies. 

Finally we have reported the results of our MD simu- 
lations of the thermal decomposition process of acetylene 
molecules on Si(lOO); in particular we observe a surpris- 
ing stability of the (dehydrogenated) C dimers and their 
tendency to diffuse below the first Si layer and to coalesce 
into more complex C fragments. 
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FIG. 1. Stable structures of acetylene adsorbed on Si(lOO), 
in the (2x1) reconstruction: B= "bridge", B'= "rotated 

bridge", P= "pedestal", P'= "rotated pedestal". For clarity 
only the four Si atoms of two dimers and four belonging to 
the second layer are shown. 

FIG. 2. Different possible structures of acetylene adsorbed 
on Si(lOO) in the B configuration, at 0.5 ML coverage: up 
atoms of Si dimers are indicated by solid circles, down atoms 
of Si dimers (and those belonging to symmetrized Si dimers) 
by open circles, and adsorbed acetylene molecules by solid 
bars. 



TABLE I. Bond length distances and vibrational frequen- 
cies of the free acetylene molecule. For comparison, besides 
the experimental values, we also report the LSDA results from 
Ref. 26 (in parenthesis). 





LSDA 


BLYP 


expt. 


C-C (A) 


1.201 (1.200) 


1.206 


1.203'' 


C-H (A) 


1.075 (1.073) 


1.066 


1.06r 


C-H stretching (cm~^) 


3304 


3337 


3374'' 


C-C stretching (crn^"'^) 


1975 


1960 


1974*= 


C-H bending (cm^^) 


739 


743 


729'' 



"■Ref. 27 
''Ref. 28 



FIG. 3. Grand canonical potential for the Si(100)-C2H2 
surface at zero temperature, as a function of the chemical 
potential McaHj ' ^'^^ ^^'^ structures (a) and (c) shown in Fig. 
2, at 0.5 ML coverage (dashed line); data are compared with 
those obtained at 1.0 ML coverage (solid line). The values are 
relative to the grand canonical potential of the clean Si(lOO) 
surface (horizontal line). 



TABLE II. Total energy per adsorbed molecule (relative 
to the B structure), E, and structural parameters for acety- 
lene on Si(lOO), at the saturation coverage of 0.5 ML. "Si" 
indicates the silicon atoms involved in bonding with acety- 
lene. Data have been computed using the gradient-corrected 
BLYP (Ref. 23) functional. Values obtained using the LSDA 
functional are instead reported in parenthesis. When slightly 
different distances are found the average value is given. See 
text for the definitions of the different structures. 

Configuration g(eV) dc-c(A) dc_si(A) dsi-si(A) 

B 0.00 (0.00) 1.36 (1.35) 1.91 (1.88) 2.37 (2.34) 

B' 0.05 (0.27) 1.37 (1.36) 1.92 (1.89) 2.43 (2.36) 

P 1.98 (1.94) 1.58 (1.54) 2.07 (2.02) 2.37 (2.33) 

P' 1.28 (0.53) 1.57 (1.53) 2.01 (1.96) 2.32 (2.28) 



TABLE III. Total energy (in rneV) per adsorbed molecule 
(relative to the (a) structure) for different structures of acety- 
lene adsorbed on Si(lOO) in the B configuration (see Fig. 2), 
at 0.5 ML coverage and zero temperature, obtained using the 
LSDA and BLYP functional, respectively. 



Configuration LSDA BLYP 

(a) 

(b) -21 89 

(c) -35 2 

(d) -22 18 
(c) -31 3 
(f) -26 17 
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